Cavity-Mediated Entanglement Generation Via Landau-Zener Interferometry 
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We demonstrate quantum control and entanglement generation using a Landau-Zener beam split- 
ter formed by coupling two transmon qubits to a superconducting cavity. Single passage through the 
cavity-mediated qubit-qubit avoided crossing provides a direct test of the Landau-Zener transition 
formula. Consecutive sweeps result in Landau-Zener-Stiickelberg interference patterns, with a visi- 
bility that can be sensitively tuned by adjusting the level velocity through both the non-adiabatic 
and adiabatic regimes. Two-qubit state tomography indicates that a Bell state can be generated 
via a single passage, with a fidelity of 78% limited by qubit relaxation. 



Avoided crossings are common in the energy level spec- 
tra of many quantum systems. In general, an off-diagonal 
matrix element that couples two states with magnitude 
A leads to an avoided crossing of magnitude 2 A in the en- 
ergy level spectrum [1]. Adjusting an external parameter 
to sweep the system through such an avoided crossing 
leads to a non-adiabatic transition probability Plz, as 
first described by Landau, Zener, Stiickelberg, and Ma- 
jorana [2]-[5]. Landau-Zener transition physics has been 
applied to atomic collisions, where the external parame- 
ter is interatomic distance, and to adiabatic rapid passage 
in nuclear magnetic resonance, where the parameter is a 
rapidly varying magnetic field O |7] . Avoided crossings 
have recently been used as effective "beam splitters" of 
quantum states in more controllable systems including 
single superconducting qubits [8HT2] and semiconductor 
double and triple quantum dots [131 E] • 

Circuit quantum electrodynamics (cQED) uses a dis- 
persive cavity-mediated interaction between two super- 
conducting qubits placed in the same microwave res- 
onator to generate an effective qubit-qubit coupling gi2 = 
/S./h [TSHlSj- The "quantum bus" architecture has en- 
abled the demonstration of two-qubit algorithms and 
three-qubit error correction [lUllO]. In this Letter, we 
use Landau-Zener physics to control states in the single- 
excitation subspace of two transmon qubits coupled to a 
superconducting resonator. Single passage through the 
avoided crossing via fast magnetic flux sweeps provides 
a direct test of the asymptotic Landau-Zener formula. 
Double passage results in tunable Landau-Zener-Stiick- 
elberg interference patterns, analogous to Mach-Zehnder 
interferometry. These experiments allow us to observe 
the full crossover from non-adiabatic to adiabatic tran- 
sitions in a single controllable system. Finally, single 
passage through the avoided crossing is used to generate 
two-qubit entanglement with a relaxation-limited fidelity 
of 78%, as verified by state tomography. Such a beam 
splitter of two-qubit states could potentially be of use in 
a quantum information processor in which tunable qubits 
utilize more than one bias point for computation. 



Our cQED device, shown in Fig. [TJa), consists of two 
aluminum transmon qubits coupled to a niobium super- 
conducting coplanar waveguide on a sapphire (AI2O3) 
substrate |T6l [iTl [211 [22] • All measurements are per- 
formed in a cryogen-free dilution refrigerator with a base 
temperature of ~ 10 mK. Each qubit is coupled to its 
own flux bias line, which we calibrate for both dc and 
ac crosstalk. Experimentally determined parameters in- 
clude a fundamental cavity mode frequency /r = 8.795 
GHz, photon decay rate /T:/27r ~ 0.5 MHz (for a cav- 
ity quality factor Q ^ 17,000), qubit-cavity couplings 
(between the first transmon transition and lowest cavity 
mode) ^i(2)/27r = 183(185) MHz, transmon charging en- 
ergies Eci{2)/h = 220(200) MHz, and transmon Joseph- 
son energies E^^^^^/h = 140(150) GHz. 

The cavity-mediated qubit-qubit coupling is first 
probed through dispersive microwave spectroscopy, 
where transmission through the cavity is measured while 
a second microwave tone is swept across the transmon 
transition frequencies [161 [23] • Transitions in the trans- 
mon qubits are detected as a dispersive shift in the cav- 
ity resonance, allowing qubit spectroscopy at frequencies 
much different from f^. Figure [TJb) reveals a coupling 
strength of 2^i2/27r = 34.8 MHz when the qubits are 
tuned into resonance at 13.9 GHz. A simple two-qubit 
gate is implemented by applying a 7r-pulse to one of the 
qubits at point I and then pulsing via the flux bias lines 
to the avoided crossing at point II for a variable amount 
of time Tp. The frequency of the resulting coherent os- 
cillations, shown in Fig. [TJc), is consistent with the en- 
ergy splitting obtained from spectroscopy. Qubit readout 
is performed using the Jaynes-Cummings nonlinearity, 
through which the strongly driven bare-cavity transmis- 
sion is highly sensitive to the state of each qubit [24| . 

To probe the Landau-Zener transition dynamics asso- 
ciated with the avoided crossing, we perform a single- 
passage experiment in which the state |10) is initialized 
at point HI and then ramped over a pulse risetime tr to 
point I, where the resulting qubit populations are mea- 
sured. A schematic illustration of the process is shown in 
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Figure 1. (color online) (a) Optical micrograph of the cQED 
device consisting of a A/2 superconducting resonator coupled 
to two transmon qubits. Qubit Q\ (Q2) is controlled with flux 
bias voltage Vi (V2) and microwave tone /i (/2). Qubit read- 
out is performed by driving the resonator with a microwave 
tone /r and measuring the transmitted homodyne voltage 
Vh- (b) Low-power spectroscopy reveals a cavity- mediated 
avoided crossing between states |10) and |01). $2 denotes the 
applied magnetic flux through the split junction loop of Q2, 
and $0 is the magnetic flux quantum. Qubit state readout 
is performed at point L (c) A time-resolved cavity- mediated 
interaction is achieved by pulsing near the avoided crossing 
(point II) for a time Tp. The blue (red) points show Vh as a 
function of Tp when the system is initialized in |10) (|01)) at 
point I. Solid lines are a guide to the eye. 

Fig.[2](a) and the specific pulse sequence is depicted in the 
inset of Fig. [2|b). The Landau-Zener formula predicts 

a non-adiabatic transition probability Plz = e~~^, 
where v is the "level velocity" of the uncoupled energies, 
V = \d[Ex - E2)/dt\ [7 . We plot the resulting popula- 
tions as a function of oc tr in Fig. [2|b). For large 
level velocities, Qi remains in its initial state and Q2 
remains in the ground state. As v decreases, Plz de- 
creases from 1, and population transfer occurs. We note 
that in the absence of relaxation, the population of Q2 
(red points) would continue to increase asymptotically 
with tr to a population of 1, corresponding to a com- 
pletely adiabatic process. Selecting tr ~ 40 ns (1/i^ ~ 13 



ns//ieV) gives a 50-50 beam splitter. Deviations from 
^LZ = 1/2 at this point are due to qubit relaxation dur- 
ing the experiment. Fitting an exponential decay to the 
Qi curve (blue points) according to the Landau-Zener 
formula (after subtracting out the independently mea- 
sured Ti^i = 120 ns Qi relaxation rate) yields a best fit 
value 2^12 /27r = 34.0 MHz, within 1 MHz of the value 
extracted from spectroscopy. Our experimental results 
agree well with numerical simulations [see solid lines in 
Fig.|2|b)] of a master equation in the Lindblad form that 
accounts for relaxation and dephasing 



ih 



(1) 

The Lindblad operator Li^q^ = y^T\^a- / describes 
relaxation of Qi to its ground state, while p2,Qi = 
\/^2,QiO'z ^ I describes pure dephasing of Qi, with the 
analogous operators defined for Q2. The weights of other 
potential Lindblad operators relevant to the dynamics 
were assumed to be small and were not included. The 
fitted rates F are consistent with coherence times ex- 
tracted from standard single-qubit Rabi and Ramsey ex- 
periments (Ti,i(2) « 120(130) ns, T*^2) ^ 150(180) ns). 
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Figure 2. (color online) (a) Fast flux bias ramps are used to 
sweep the system through the cavity- mediated avoided cross- 
ing, with dynamics governed by the Landau-Zener transition 
formula. Passage through the avoided crossing is equiva- 
lent to a beam splitter. If the system is initialized in state 
1 10) and then swept through the avoided crossing, the prob- 
ability to remain in |10) is given by Plz- (b) The cavity- 
mediated beam splitter is characterized by sweeping through 
the avoided crossing from point III to point I for a range of 
ramp times tr oc 1/z/. Qi {Q2) populations are plotted as a 
function of l/z^ in blue (red). 
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Figure 3. (color online) (a) A sweep back and forth across the avoided crossing is equivalent to a Mach-Zehnder interferometer, 
with level velocity-dependent transmission ratios that are controlled by the ramp time tr of the flux bias pulse. The red and 
blue curves depict how the eigenenergies evolve during the experiment. The labels I and II on the left refer to the labeled flux 
bias points in Fig. [lib), (b) P|io) plotted as a function of flux bias pulse width tp and flux bias pulse detuning A$2 for ramp 
times Tr ~ 1.5 ns (left) and tr = 18 ns (right). A pulse to the avoided crossing (point II) corresponds to |A$2| ~ 12 m$o- 
Insets: Simulated Landau-Zener-Stiickelberg interference patterns, (c) P|io) as a function of rp for tr 1.5 ns and tr = 18 ns, 
as extracted from the data sets in panel (b) (dashed lines), (d) Oscillation visibility V plotted as a function of l/u. The green 
solid line is the result of the double beam splitter (DBS) formula (jsj after plugging in previously extracted parameters and 
neglecting relaxation. The solid red line is the result of multiplying the green curve by an exponential decay (with previously 
measured time constant Ti,i ^ ^^1,2) with respect to l/u oc tr, with offset and scaling based on measured ramp and buffer 
times (see Supplemental Material), and agrees closely with a curve obtained using the numerical master equation simulation. 



The Landau- Zener formula describes transition proba- 
bilities, but not the more fundamental transition ampli- 
tudes that give rise to many interesting quantum interfer- 
ence phenomena. Phase is an important relation between 
two interacting states, in particular states that "recom- 
bine" after a beam splitter event. If the avoided cross- 
ing is doubly traversed (across and back), the relative 
phase accumulated between the two consecutive cross- 
ings will lead to Stiickelburg oscillations [4 . In conven- 
tional Landau- Zener problems such as atomic collisions, 
phase accumulation occurs so rapidly that interference is 
washed out by even a small amount of decoherence [4 . 
However, in artificial atoms this phase can be observed 
and precisely tuned [9} [TO l IT2 l 113] . 

In Fig. |3] we observe Stiickelburg oscillations by per- 
forming a "double beam splitter" experiment that is anal- 
ogous to Mach-Zehnder inter ferometry. We are able to 
sensitively tune the visibility of the resulting Stiickel- 
berg oscillations, in agreement with theoretical predic- 
tions. Schematics of a double passage experiment in 
the non-adiabatic (Plz ^ 1) and perfect beam splitter 
{Plz ^ 0.5) regimes are shown in Fig. [3]^ a). State |10) 



is prepared near point I and linear flux bias ramps are 
applied to sweep the system through the avoided cross- 
ing (point II) and back, with symmetric ramp time tr 
and pulse duration rp. The resulting experimental data 
are shown in Fig. [sj^b), again with the left panel having 
Plz ^ 1 (tr ~ 1.5 ns, leading to ~ 1 /ieV/ns and 
Plz ^ 0.95 at the dashed line) and right panel having 
Plz ^ 0.5 (tr = 18 ns, leading to v ^ 0.08 /ieV/ns 
and Plz ^ 0.56 at the dashed line). The resulting in- 
terference patterns agree well with the numerical master 
equation simulation ([T]), as shown in the insets of Fig. 
[3jb). Horizontal cuts through these data are shown in 
Fig. [sjc), indicating that the visibility of the oscillations 
is a sensitive function of the level velocity. Note that 
these oscillations occur in their entirety below P|io) = 1 
because of qubit relaxation during ramp and buffer times. 

The level velocity-dependence of the oscillation vis- 
ibility can be understood by examining the successive 
transition amplitudes during the double passage experi- 
ment. Neglecting relaxation and dephasing, each passage 
through the beam splitter performs an effective unitary 
operation on the incoming state according to the transfer 
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matrix [1011171 [H] 



U 




(2) 



where (^s = — 7r/2 + ^s, 0s = 7r/4 + Arg [r(l — i5)] + 
(5(ln5 — 1) is the Stokes phase, and V is the gamma func- 
tion. The Stokes phase is a function of the adiabaticity 
parameter 5 = [IH]- Inserting a free evolution pe- 

riod between the two beam sphtter operations results in 
a return probability 



|10) 



1-2Plz(1-Plz) 1 + cos((/)-2^s) , (3) 



where (j) is the dynamical phase accumulated on the far 
side of the avoided crossing. Therefore, in the absence 
of decoherence, the visibility of the Landau-Zener-Stiick- 
elberg interference fringes is V = 4Plz(1 — Plz)- We 
explicitly test this dependence by measuring coherent os- 
cillations as a function of rp at a fixed |A<I>2| ~ 16 m<l>o 
for different beam splitter ramp times tr. The resulting 
oscillation visibilities, defined as the initial amplitudes of 
fits to the decaying oscillations (see Supplemental Mate- 
rial), are plotted as a function of 1/v in Fig. jsjd). Our 
results agree well with theory when relaxation is taken 
into account and yield a maximum oscillation visibility 
when 1/v ^ 12 ns//ieV, consistent with the single-pass 
beam splitter calibration in Fig. [2|b). 

Finally, we demonstrate entanglement generation 
through a single passage in the 50-50 beam splitter limit. 
After initializing the system in state |10) at point III, 
we sweep across the avoided crossing with a ramp rate 
corresponding to Plz =0.5. This time, a set of tomog- 
raphy pre-rotations /7k G SU (2) SU (2) are performed 
at point I directly before measurement. Figure |4ja) il- 
lustrates the ideal process in the absence of relaxation 
and dephasing. The initial state = |10) is inci- 

dent on the beam splitter, yielding an entangled state 
IV^out) = 1/V2(|10)-|-|01)). The results from state to- 
mography are shown in Fig. [4[b ) and agree well with the 
master equation simulation ^TJ that takes into account 
relaxation and dephasing (red dashed lines). We obtain 
a state fidelity with respect to the target Bell state of 
F = (V^bIpIV^b) = 78% and entanglement of formation 
Ei = 63% [29], both limited by relaxation. 

In conclusion, we have utilized cavity-mediated 
Landau-Zener physics to achieve coherent control of a 
two-qubit transmon system. We have verified the general 
Landau-Zener formula and the dependence of double- 
passage Stiickelberg oscillations on the level velocity. 
A simple entanglement generation protocol is demon- 
strated, where a single passage through the avoided cross- 
ing in the 50-50 beam splitter limit yields a Bell state 
with relaxation-limited F = 78%. The combination of 
single qubit rotations and Landau-Zener two-qubit in- 
teractions may be used to implement more complicated 
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Figure 4. (color online) (a) In the 50-50 beam split- 
ter limit, a single pass through the avoided crossing con- 
verts the input state \ipin) = |10) to the Bell state iV^out) = 
1/^(|10) + |01)) = I^b). (b) An entangled state is gener- 
ated by preparing the system in state |10) and then sweep- 
ing through the avoided crossing with a ramp time corre- 
sponding to l/iy — 13 ns//ieV and Plz = 0.5. We per- 
form state tomography to extract the density matrix of the 
two-qubit system via maximum-likelihood estimation [19] , 
assuming a joint measurement operator of the form M — 
Pull + hiZI + PizIZ + hzZZ 30 . The density matrix 
shown is after an azimuthal Z rotation on Q2 chosen to max- 
imize the off-diagonal terms in Re [p] . Imaginary parts (not 
shown) are all small compared to the real components (less 
than 0.04). See Supplemental Material for further experimen- 
tal details. 



quantum algorithms in systems with longer coherence 
times. 
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